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Abstract 

We investigate the behavior of the Shanghai Stock Exchange Composite (SSEC) index 
for the period from 1990:12 to 2007:06 using an unconstrained two-regime threshold au- 
toregressive (TAR) model with an unit root developed by Caner and Hansen. The method 
allows us to simultaneously consider non-stationarity and nonlinearity in financial time se- 
ries. Our finding indicates that the Shanghai stock market exhibits nonlinear behavior with 
two regimes and has unit roots in both regimes. The important implications of the threshold 
effect in stock markets are also discussed. 
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1 Introduction 



In the past three decades since the Third Plenary Session of the 1 1th Central Com- 
mittee of the Communist Party of China in December 1978, China has paved a 
gradual transition from a centrally planned economy to a market economy and the 
economy has experienced unprecedented growth. During this period, one of the 
most important developments has been the reopening and operation of the Chinese 
stock market. Before the foundation of People's Republic of China, the Shang- 
hai Stock Exchange was the third largest worldwide (after New York and Lon- 
don Stock Exchanges) and had remarkable influence on other world-class financial 
markets fP|. After 1949, China implemented policies of a socialist planned econ- 
omy and the government controlled entirely all investment channels. In 1981, the 
central government began to issue treasury bonds to raise capital to cover its fi- 
nancial deficit, which reopened China's securities markets. The first market for 
government- approved securities was founded in Shanghai on November 26, 1990 
and started operating on December 19 of the same year under the name of Shanghai 
Stock Exchange (SHSE). Shortly after, the Shenzhen Stock Exchange (SZSE) was 
established on December 1, 1990 and started its operations on July 3, 1991. The 
size of the Chinese stock market has increased remarkably |l2l. There are increas- 
ing interests in the academic studies of the Chinese stock market. 

It is well-known that most econometric models are constructed based on the as- 
sumption that the variables are stationary. It is thus very important to perform unit 
root test for stationarity. Empirical studies on the US markets find mixed results 
concerning the presence of a unit root in the behavior of stock indexes [|3l. The 
situation seems alike in the Chinese stock market. Xie, Gao and Ma applied the 
Phillips-Perron (PP) test to the weekly data of the Shanghai Stock Exchange Com- 
posite (SSEC) index from 12/21/1990 to 03/02/2001 and found that the nuU hypoth- 
esis that the logarithm of the index has a unit root cannot be rejected at the signifi- 
cance level of 5% BUl. Cheng, Wu and Zhou tested the unit root property in the daily 
data of SSEC from 01/02/1998 to 12/31/2001 using the augmented Dickey-Fuller 
(ADF) approach and reached a similar conclusion |[5l. In contrast, Dai, Yang and 
Zhang investigated the daily data of the SSEC index from 12/19/1990 to 06/1 8/2004 
with the ADF approach and found that the unit root null is rejected at the 5% level 
of significance Q. 

On the other hand, numerous evidence shows that there exists threshold nonlinear- 
ity in the behavior of stocks II7I8II . It is thus helpful to distinguish non- stationarity 
from nonlinearity in the stock market behavior. This task can be done by adopting 
the threshold unit root test developed by Caner and Hansen [9]. Briefly speaking, 
Caner and Hansen use threshold autoregressive (TAR) model to test for a threshold 
effect and then perform unit root tests on both regimes if exist. Recently, Narayan 
applied this method to the stock prices of Australia (ASX All Ordinaries, monthly 
data from 1960:01 to 2003:04) and New Zealand (NZSE Capital Index, monthly 
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data from 1967:01 to 2003:04) ITOl, and the US stock price index (NYSE common 
stocks, monthly data from 1964:06 to 2003:04) [31. The main finding is that the 
stock prices in the three markets are generated by nonlinear processes and can be 
characterized by unit root processes. 

This work attempts to add to the existing literature by testing for nonlinearity and 
unit root property of the Chinese stock price index SSEC, which is monthly over 
the period from 1990: 12 to 2007:06. We find that the monthly SSEC index behaves 
nonlinearly with a unit root in both regimes. The paper is organized as follows. 
Section [2] reviews briefly the procedure of the threshold unit root test of Caner and 
Hansen Section [3] presents the empirical results. And Section |4] gives some 
conclusive remarks. 



2 Econometric methodology 

In this section, we describe briefly the econometric methodology of the threshold 
unit root tests proposed by Caner and Hansen A rigorous presentation with 
assumptions, theorems and proofs can be found in their seminal paper [9|. See also 
references II11I12II for a tutorial example. 

2. 1 Model specification and calibration 

Following the work of Caner and Hansen dU, we adopt a two-regime threshold 
autoregression (TAR) model with an order of k. The mathematical expression of 
the TAR(/!;) model reads 

Ayt = e,Xt_J{Zt_, < A) + ^2Xt-i/(Z,_i ^ A) + Q (1) 

with 

xt^i = ivt-i, 1, Ayt^i, ■ • ■ , Ayt-k)' , (2) 

where y is the logarithm of the SSEC index for t = 1, 2, • ■ • , T, is an i.i.d. 
error, /(expression) is the indicator function that equals to 1 if the expression in 
the parentheses is true and otherwise, Zt = yt — yt~m for some m. ^ 1 is the 
threshold variable, and A; ^ 1 is the autoregressive order. The variable Zt has clear 
financial meaning acting as return at the time horizon of m months. The threshold 
parameter A is unknown and represents the level of the variable yt that triggers a 
"regime change", if any. The components of 9i and 62 can be partitioned as follows: 

\0i = (pi,/?i,ai) 
I O2 = {92,^2,^2) 
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where pi and p2 are slope coefficients on yt-i. Pi and (32 are scalar intercepts, and 
ai and Q!2 are 1 x A; vectors containing the slope coefficients on dynamics regressors 
(A2/t„i, . . . , Ayt^k) in the two regimes. 



In order to calibrate model ©, the concentrated least squares approach is usually 
utilized. The regression procedure is carried out for each value of m. The value of 
A is taken from a compact interval [Ai, A2] in which Ai and A2 are determined by 
the following constraints 

' Pr(Z, ^ Ai) = TTi 

Pr(Zj ^ A2) = TT2 ' 



where < tti < 7r2 < 1 and tti + 712 = 1. In this work, we impose tti = 0.15 lUTll . 
For each AG [Ai, A2], the parameters p's, /9's and a's are estimated by minimizing 
the objective function 

T 

Qi\,m) = Y.etiX,mf . (5) 

t=i 

Let et{X, m) represents the residual from the ordinary least squares for given A and 
m. Then the least squares estimate A of the threshold parameter is given by 

A = min Q(X, m) . (6) 

Ae[Ai,A2]^' 

Note that A and other estimates of parameters are dependent of m. 



2.2 Test for threshold effect 



In model a question of particular interest is whether or not there is a threshold 
effect. The threshold effect disappears under the null hypothesis 

Ho -.9^ = 92, (7) 

which is tested using a standard heteroskedastic-consistent Wald test |l9l. The Wald 
statistic is 

W = W{X) = sup W{X) . (8) 

AG[Ai,A2] 

If the null hypothesis can not be rejected, there is no threshold effect, in which case 
the two vectors of coefficients are identical between the two regimes (^1 = 6*2). 
Caner and Hansen also show that sup;^^^;^^ ^2] W{X) has a non-standard asymptotic 
null distribution and propose a bootstrap method to compute the asymptotic critical 
values and p- values BUl. 
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2. 3 Tests for threshold unit root 



When there are two regimes delimited by a threshold, we have two parameters pi 
and p2 controlling the stationarity of the process yt. The null hypothesis is 

Ho:pi = p2 = 0. (9) 

When the null hypothesis Hq holds, the process yt has a unit root and model ([T]) can 
be expressed in terms of the stationary difference Ayt. An alternative hypothesis to 
the null Hq is 

Hi: pi<0 and P2 < . (10) 

When Hi holds, the process yt is stationary and ergodic in both regimes [|9]|. An- 
other alternative deals with a partial unit root, which is expressed as follows 

Pi < and p2 = 0, or 
H2 ■ \ (11) 
Pi = and P2 < 0. 

When H2 holds, the process yt have a unit root in one regime and is stationary in 
the other showing mean reversion behavior. 

To test the null hypothesis Hq against its two alternatives Hi and H2, there are two 
Wald tests that apply. The statistic of one-sided Wald test against the unrestricted 
alternative pi 7^ or p2 7^ is 

Ri=tllipi<0)+tll{p2<0), (12) 

while that of the two-sided Wald test against pi < or p2 < is 

R2 = tl + tl, (13) 

where ti and ^2 are the t ratios for pi and p2. We note that this term t should not be 
confused with the time t in Eq. ([T]). In order to further discriminate the stationarity 
in the two regimes, we can examine the negative of the t statistics —ti and — ^2- 



3 Application to the SSEC index 



3.1 The data set 



We analyze the whole time series of the SSEC index. The components of the SSEC 
index consist of all stocks listed on the Shanghai Stock Exchange, including both 
A shares and B shares. It is thus an overall index reflecting the price fluctuations 
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of the overall Shanghai stock market. The index was officially released since July 
15, 1991, tracing back to December of 1990. Monthly data over the period from 
1990:12 to 2007:06 are utilized for analysis. Specifically, we retrieve the closing 
prices of the last trading days of all months and calculate the logarithms, which 
give the time series yt defined in the preceding section. 



3.2 Unit wot tests 



As a first step we perform conventional unit root tests of the monthly SSEC in- 
dex without taking into account possible nonlinearity. The augmented Dickey- 
Fuller (ADF) [m, Phillips-Perron (PP) flT], and Kwiatkowski-Phillips-Schmidt- 
Shin (KPSS) [[T5l tests are adopted. We test for unit root in both the logarithm 
of SSEC yt and its first-order difference Ayt. For the ADF and PP tests, the null 
hypothesis is that yt (resp. Ayt) has a unit root, which utilizes the t-statistic. In 
contrast, the null of the KPSS method is the stationarity of the variable and uses the 
LM-statistic. The results are presented in Table [B All tests indicate that the SSEC 
index is a unit root process while its first-order difference is stationary. 

[Table 1 about here.] 



3.3 Threshold effect 

We now use the Wald test to examine whether we can reject the linear autoregres- 
sive model in favor of a threshold model. In our model, we adopt that k = 12. In 
Table [21 we report the results of the Wald test. Also listed are the bootstrap critical 
values at three conventional levels 10%, 5%, and 1% and the bootstrap p-values for 
threshold variables of the form Zt = yt — yt~m for different delay parameters m 
ranging from 1 to 12. The bootstrapping is carried out with 1000 and 2000 replica- 
tions. The results are qualitatively the same for both cases so that we report below 
the results with 2000 replications. For all m, the null hypothesis 9i = 62 of lin- 
earity is rejected at the significance level of 1%. In other words, the presence of a 
threshold effect in the monthly SSEC index is statistically significant with a 99% 
confidence level. According to these results, the linear AR model can be rejected in 
favor of the TAR model. 

[Table 2 about here.] 

The optimal value of decay m can be determined exogenously, which maximizes 
the value of |[9l. According to Tabled the Wald statistic is maximized (W = 
127.1) when m = 4. Hence, we take m = 4 as the optimal decay parameter, which 
results in a preferred TAR model. Accordingly, the point estimate A of the threshold 
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is determined to be —0.1418. For the preferred specification with m = 4, we report 
in Table |3] the least squares parameter estimates 9i and 62 with standard errors. 

[Table 3 about here.] 

The TAR model identifies two regimes depending on whether the variable Zt = 
yt — yt-4 lies above or below the threshold A = —0.1418. The first regime is when 
Zt < —0.1418, which occurs when the SSEC index has fallen cumulatively more 
than 14.18% in the last four months. About 16.4% of the observations fall into this 
first regime. The second regime is for Zt ^ —0.1418, which constitutes all those 
observations that occur when the m-month price variation is no less than —0.1418. 
Approximately 83.6% of the observations belong to the second regime. Figure [T] 
shows the estimated division of the SSEC index into two regimes. 

[Fig. 1 about here.] 

3.4 Threshold unit root tests 



We examine the unit root properties of the SSEC index that possesses significant 
threshold effect. We first compute the one-sided and two-sided threshold unit root 
test statistics Ri and R2 together with the bootstrap critical values at three signif- 
icance levels 10%, 5% and 1% and p- values for each delay parameter m, ranging 
from 1 to 12. The critical Wald-statistics at three significance levels as well as the 
p-values are calculated according to a bootstrap approach with 2000 replications. 
The results are reported in Table HI The one-sided Wald tests in the left panel of 
Table m show that the statistic Ri is less than Wio%. The situation is similar for the 
two-sided Wald tests presented in the right panel of Table |4l In summary, for all 
m, both i?i and R2 are less than the critical value at the 10% level of significance. 
These results suggest that the null hypothesis of the presence of a unit root in the 
monthly SSEC index cannot be rejected at the 10% level of significance. 

[Table 4 about here.] 

Although both tests Ri and R2 cannot reject the unit root hypothesis, they are not 
able to discriminate between the full unit root case in both regimes and the partial 
unit root case in one regime. We thus test the partial unit root in the monthly SSEC 
index by calculating the individual t statistics, and ^2- The results are reported 
in Table |5l The critical Wald-statistics at three significance levels as well as the 
p-values are calculated according to a bootstrap approach with 2000 replications. 
We find that, for all m, both ti and ^2 are less than the critical value at the 10% 
significance level. Hence, we are again unable to reject the unit root null hypothesis 
in both regimes of the monthly SSEC index. 

[Table 5 about here.] 
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It is noteworthy that same conclusions are reached when we use in the above sta- 
tistical test the asymptotic valued tabulated by Caner and Hansen BU. 



4 Concluding remarks 

In summary, we have adopted the econometric approach of threshold autoregres- 
sion (TAR) with a unit root developed by Caner and Hansen ||9l to analyze the 
monthly data of the Shanghai Stock Exchange Composite index. The SSEC index 
is found to have a threshold effect of A = —0.1418 with strong evidence. In ad- 
dition, both regimes with the index variation below or above the threshold have 
significant unit roots, so does the whole time series. Our results indicate that the 
Shanghai stock market exhibits nonlinear behaviors with a unit root. 

An important question arises asking what we can learn further from the fact that 
the stock market is nonlinear with a threshold. The presence of a threshold A = 
—0.1418 means that the market behaves differently when it falls more than 14.18% 
in four months. This threshold effect has direct connection with the concept of large 
drawdowns in the sense of coarse graining in time for the former and price variation 
for the latter, which are usually outliers I116I17L By scanning different time scales, 
one might be able to provide evidence for such a connection. 

Another closely relevant issue is the definition of crashes. A consensus is still 
lack. A quite feasible and unambiguous definition is based on large drawdowns 
[1181191 . A systematic investigation shows that more than 50% of the crashes iden- 
tified are endogenously triggered that have evident preceding log-periodic power- 
law (LPPL) patterns [12012 111 . An alternative option is to seek for large price drops 
within different time windows [|22]| . These two methods identify partially overlap- 
ping examples of crashes. It is thus interesting to explore the possibility of having 
a new definition of crashes by developing a multiscale TAR approach. The idea 
of multiscale analysis is also related to the LPPL investigation of large financial 
variations [|23l . It is however beyond the scope of the current work. 

Similar to the parity pair of drawdown and drawup, it is natural to think of the pres- 
ence of a positive threshold. This needs a three-regime threshold autoregression 
method with two thresholds. If there do exist two thresholds (positive and negative) 
in the stock market behavior, one can expect that the two-regime threshold autore- 
gression method will result in a negative threshold in some cases as in the monthly 
SSEC index and a positive threshold in other cases. This calls for further studies. 
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Fig. 1. Monthly data of the SSEC index classified by threshold regime. The first regime 
("Regime 1" in the legend) constitutes observations with the cumulative price drop greater 
than 14.18% in the last four months. 
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Table 1 

Unit root tests of the monthly SSEC index without thi^eshold effect. The augmented Dickey- 
Fuller (ADF), Phillips-Perron (PP) and Kwiatkowski-Phillips-Schmidt-Shin (KPSS) tests 
are adopted. The left panel tests for unit root in the logarithm of SSEC yt, while the right 
panel tests for unit root in the first-order difference. Vx% is the critical value at the x% 
significance level. 

log SSEC, j/, A log SSEC, Ay^ 

method statistic Vio% ^^5% Vi% j»-value statistic 1^10% Vi% p-value 

ADF -0.76 -1.62 -1.94 -2.58 0.877 -2.56 -1.62 -1.94 -2.58 0.011 
PP 1.55 -1.62 -1.94 -2.58 0.970 -14.6 -1.62 -1.94 -2.58 0.000 
KPSS 1.28 0.35 0.46 0.73 0.000 0.16 0.35 0.46 0.73 >0.1 
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Table 2 

Wald tests for a threshold effect in the monthly SSEC index for different lags m. The 
second row gives the Walt-statistics W for different m. The third to fifth rows show the 
critical Wald-statistics at three significance levels according to a bootstrap approach with 
2000 replications. The last row presents the bootstrap p- values. The optimal delay is m = 4, 
highlighted in bold face. 



m 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


W 


84.0 


100.2 


120.4 


127.1 


88.4 


63.7 


68.7 


102.2 


113.0 


102.9 


87.2 


110.2 




44.1 


42.7 


41.8 


43.0 


41.8 


42.0 


40.1 


41.5 


41.2 


41.1 


39.3 


39.9 




54.7 


50.9 


51.4 


52.2 


53.9 


52.6 


51.0 


50.4 


51.5 


51.6 


48.6 


51.8 




75.5 


72.6 


76.5 


72.5 


75.6 


81.1 


76.7 


75.5 


77.2 


77.2 


74.8 


74.7 


p-value 


0.4% 


0.1% 


0.1% 


0.1% 


0.4% 


2.4% 


1.8% 


0.1% 


0.1% 


0.1% 


0.4% 


0.1% 
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Table 3 

Least-squares estimates of parameters from the unconstrained threshold model with an 
optimal decay m = 4. The threshold estimate is A = —0.1418. 



Regressors 




< A 


Zf 1 


> A 


Estimate 


S.E. 


Estimate 


S.E. 


i/i— i 


0.0091 


0560 


-0.0109 


0.0197 


Tntercent 


-0.1272 


0.4392 


0.0793 


0.1414 




-1.9509 


0.3440 


0.1530 


0.0705 




0.0504 


0.5108 


-0.0495 


0.0671 




1.0005 


0.4000 


0.0978 


0.0671 


A?/* /I 
'-^yt—i 


-0.6906 


0.3871 


0.0379 


0.0634 


'—^yt—o 


0.5264 


0.1916 


0.0852 


0830 


A?/+ fi 
'—^yt—b 


-0.1085 


0.1564 


0.0415 


0.0837 




0.5258 


0.1667 


-0.0564 


0.0983 




-0.6614 


0.2292 


-0.0318 


0.0658 


^yt~9 


1.4061 


0.2966 


0.2429 


0.0626 




1.1534 


0.3654 


-0.0484 


0.0656 


^yt-ii 


-0.5198 


0.3796 


0.2115 


0.0634 




0.8813 


0.2942 


-0.0290 


0.0668 
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Table 4 

One-sided (left panel) and two-sided (right panel) Wald tests for threshold unit roots in the 
monthly SSEC index for different lags m. The optimal delay is m = 4, highlighted in bold 
face. 



One-sided Wald test, Ri Two-sided Wald test, R2 



m 


Ri 


Wio% 




w,% 


p-value 


R2 


Wio% 


W5% 




jo- value 


1 


0.9 


11.2 


14.6 


23.9 


0.838 


4.2 


11.5 


15.0 


23.9 


0.511 


2 


1.1 


11.2 


14.1 


22.8 


0.819 


8.3 


11.6 


14.5 


23.0 


0.216 


3 


4.5 


11.9 


15.6 


27.2 


0.446 


4.5 


12.6 


16.4 


27.8 


0.496 


4 


0.3 


11.7 


15.4 


25.5 


0.922 


0.3 


12.1 


15.6 


25.7 


0.975 


5 


0.0 


12.2 


16.6 


31.5 


0.973 


0.0 


12.5 


17.5 


31.5 


0.997 


6 


1.0 


12.0 


16.4 


28.0 


0.845 


1.0 


12.5 


16.5 


28.9 


0.901 


7 


3.4 


12.5 


16.4 


30.2 


0.587 


3.4 


12.9 


16.7 


30.5 


0.647 


8 


0.1 


13.1 


17.2 


28.2 


0.954 


0.2 


13.5 


17.6 


28.6 


0.978 


9 


1.4 


12.9 


17.4 


33.3 


0.805 


1.4 


13.5 


18.0 


33.9 


0.863 


10 


1.0 


13.0 


17.6 


32.0 


0.857 


1.7 


13.5 


18.0 


32.5 


0.845 


11 


0.4 


13.8 


18.1 


35.4 


0.919 


0.7 


14.4 


19.1 


37.4 


0.935 


12 


0.4 


14.7 


19.7 


42.6 


0.911 


1.9 


15.0 


19.9 


42.6 


0.815 
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Table 5 

Wald tests for thieshold unit roots in the two regimes of the monthly SSEC index for dif- 
ferent lags m. The optimal delay is m = 4, highlighted in bold face. 



TTl 


i-siai 








^Vk I/O 111 £!i 

j9-vdiue 


cf of 

t-Sldl 








p- value 


1 
1 


1 c 

-1.0 


Z.D 


-5.U 


A 'J 

4. J 


u.y / / 


u.y 


z. / 


'J 

J.z 


A 
4.Z 


U.J 


1 


-1. 1 


Z.D 


J.U 


A 1 

4.1 




l.U 


Z. / 


J.z 


A 

4.Z 


U.jzj 


i 


z.U 


Z.D 




A A 

4.4 


U.zU/ 


U. / 


Q 
Z.O 


J. 4 


A f. 
4.D 


U.DZO 


A 

4 


-U.Z 


Z.3 


1 A 


/I 1 
4.1 


U.oZ/ 


U.o 


Z.I 


1 1 


4. / 


U.O /J 


c 

J 


-U.Z 


Z.D 


'J 1 


A ^ 

4. J 


U.OZD 


-U.U 


Q 
Z.O 


'J < 
i.j 


J. 1 


U.olo 


6 


1.0 


2.6 


3.2 


4.2 


0.522 


0.3 


2.9 


3.5 


4.9 


0.737 


7 


1.8 


2.6 


3.1 


4.1 


0.277 


0.0 


2.9 


3.5 


5.1 


0.812 


8 


-0.3 


2.7 


3.2 


4.4 


0.871 


0.3 


2.9 


3.6 


4.9 


0.724 


9 


0.5 


2.7 


3.2 


4.3 


0.708 


1.1 


3.0 


3.6 


5.6 


0.517 


10 


-0.8 


2.6 


3.1 


4.2 


0.921 


1.0 


3.0 


3.6 


5.4 


0.554 


11 


-0.6 


2.8 


3.3 


4.5 


0.893 


0.6 


3.0 


3.7 


5.6 


0.658 


12 


-1.2 


2.8 


3.4 


4.6 


0.952 


0.7 


3.1 


3.8 


6.1 


0.653 
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